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Abstract 

The spin-1 XY chain in a transverse field is studied using finite-size scaling. The ground 
state phase diagram displays a paramagnetic, an ordered ferromagnetic and an ordered os- 
cillatory phase. The paramagnetic-ferromagnetic transition line belongs to the universality 
class of the 2D Ising model. Along this line, universality is confirmed for the finite-size 
scaling functions of several correlation lengths and for the conformal operator content. 
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In modern theories of (equilibrium) critical phenomena, the notions of scaling and universality 
play a central role. These notions are particularly useful when applied to finite systems using 
finite-size scaling techniques, see ^ for an extensive review. In this work, we study the effects of 
varying the spin quantum number on the thermodynamics of the well-known XY quantum chain 
in a transverse field. For spin |, this model is exactly integrable in terms of free fermions and 
many of its properties are well studied, see §]. Besides being of interest in its own right (i.e. for 
the influence of the quantum effects on the order parameter profile [^), this quantum Hamiltonian 
also arises from the master equation description of several non- equilibrium statistical systems, see 
P]. Here we consider the spin-1 variant of this model, with the Hamiltonian 
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where h is the transverse field, 77 measures the spin anisotropy, C is a normalization constant and 
N is the system size. We use periodic boundary conditions. Finally, the S^'^'^ are spin-1 matrices 



5" = (2) 






(Spin-1 Ising models were recently proposed to describe the adsorption of CO on graphite, see 
0.) We are interested in the ground state energy Eq, which plays the role of the equilibrium 
free energy (for reviews see 0, §|) and in the correlation lengths related to the exponential 
decay of two-point correlation functions, given by the energy gaps ^^^^ = Ei — Eq. We calculate 
the low-lying spectrum of H for finite (up to = 14) using the Lanczos algorithm and then 
extrapolate towards N 00, see for details. The quantum Hamiltonian H commutes with the 
charge operator Q, the parity operator P and the translation operator T defined by 

N 

g = n {ns-r - 1) , ps^^y'-p^ = 5^ti-n , Ts^'y'-T^ = s:tf (3) 

i=l 

Eigenstates of H are thus characterized by the eigenvalues of Q, P and T, which serves to block- 
diagonalize H. 

Our first task is to determine the phase diagram, shown in Fig. 1. We recognize three distinct 
phases. The first transition, between the paramagnetic phase P and the ferromagnetic phase F, 
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Table 1: Critical points hc{i]) and conformal normalization ({r]) for the spin-1 XY model along 
the Ising line. The numbers in brackets give the estimated uncertainty in the last digit. 



is found from conventional finite-size scaling and will be shown below to be in the 2D Ising 
universality class. Close to a conventional critical point of second order, the following finite-size 
scaling form for the inverse correlation lengths is expected ||^, [T| 

= N-'S, {CNy{h - K)) (4) 

where he = hc{r]) is the critical point, y = 2 — a. critical exponent, C is a non- universal metric 
factor and Si is a universal scaling function. In particular, from 2D conformal invariance, it follows 
that Si{0) = 2Txxi [0], where Xi is a universal critical exponent. Now, the critical point h^ can 
be found from phenomenological renormalization 0. The results, extrapolated to — oo, are 
displayed in table |I[ For = 1 we find agreement with the earlier result []1T| he — 1.3259. 



The second transition occurs between the ferromagnetic phase F and a new 'oscillatory' 
phase O. This transition is well known for the spin-| case [0 and occurs along the line h = hoirf) 
where0 

+ Ki'nf = 1 (5) 

For spin- 1 it is known that while in the F phase the connected spin- spin correlation function 
< S^S'q >e decays monotonously with R, the oscillatory phase is characterized by a new wave 
vector K which modulates the spin-spin correlator 

< S^j^S'^ exp(-2i?/0 cos{KR) (6) 

Furthermore, in the oscillatory phase there are level crossings in the ground state energy which 
occur already for finite values of the number of sites |l^, ^. It was shown in |jl2| that the 
location hk{N) of the k-th level crossing satisfies a finite-size scaling law 

he-hk{N)r^N-'/'' (7) 

^ Along this line H can also be obtained from the master equation of certain ID stochastic systems ||]. 
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Table 2: Extrapolated finite-size estimates for the critical point ho and the exponent v for the 
spin-1 XY model (|I|) as determined from the second ground state level crossing h2{N) of eq. (|^. 



where the exponent u describes the scaling of the wave vector K ~ {ho — hy in the vicinity of the 
0/F transition line (for h < ho )■ For spin ^ = |, it is known that u = 1/2 [|12[. 

We now ask whether a similar transition occurs for larger values of S. Indeed, it is known that 
for arbitrary spin S and periodic boundary conditions, the ground state energy of H is doubly 
degenerate aX h = ho{ri) |13|. For spin S" = 1, we have checked numerically that the first ground 
state level crossing hi{N) always occurs at h = ho{T]) for all finite N. In addition, the second 
crossing h2{N) converges towards ho{ri), as apparent from the extrapolated data in table |[ The 
exponent found from eq. (|^) is consistent with u ~ 1/2, independently of rj and in agreement 
with the exact result for spin-|. This supports universality along the F/0 transition line. In fact, 
having confirmed the same finite-size scaling behaviour of the level crossings in the ground state 
energy for both spin S" = | and S" = 1, we expect the features of the oscillatory phase known 
P, |12| from 5 = I to be present for S* = 1 as well. 

From now on, we concentrate on the P/F transition line. We expect this transition to be in the 
2D Ising universality class, if rj ^ 0. To see this, we compare the low-lying excitation spectrum of 
H with the prediction of conformal invariance |T^, |^, following the steps outlined for the spin-| 
case in H, p. 135]. Conformal field theory states that, after subtraction of a purely extensive 
term, H can be written in the form 

where c is the central charge and Lq, Lq are generators of the Virasoro algebra which acts as 
a dynamical symmetry for H. As a consequence, eigenstates can be grouped into 'conformal 
towers', each represented by exactly one primary operator with conformal weights (A, A). The 
scaling dimension of the corresponding eigenstate is a; = A-|- A. The scaled energies and momenta 
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take the form 

^a,a(/, /) = lim (e^A^, I) -Eo)-^ = {A + I) + iA + I) (9) 

^a,a(/, /) = lim Pa,a(/, /) ■ ^ = (A + /) - (A + 7) (10) 

with I, I integer. -Ea^A; -Pa.A; respectively, are the eigenvalues of H,i\nT and Eq is the ground 
state energy. However, the application of these relations requires that the scaled energies S and 
momenta V are measured in the same units, thus fixing the normalization ( of H accordingly. We 
find ( by demanding that £^o,o(2, 0) = 2 throughout |]14|. The results for ( are given in table [|. 



Next, we determine the central charge. For = 1, we find c = 0.49999(1), close to the expected 
c = I for the 2D Ising universahty class. We did not compute c explicitly for other values of r], 
but expect c to be 77- independent. In order to check the complete operator content, we give the 
extrapolated values of the scaled energies in the charge sectors Q = and Q = 1 in tables |^ and |. 
When comparing these spectra to the expected operator content of the 2D Ising model |]14|, |^, 
namely for the Q = sector the conformal towers generated by the primary operators (0, 0) and 



|,|) (which correspond to the vacuum 1 and the energy density e) and for the Q = 1 sector the 

1 1 



conformal tower generated by (j^ij^) (which corresponds to the order parameter density cr), we 
find complete agreement. In particular, we read off the scaling dimensions x^- = | and = 1 
which determine the bulk critical exponents. 

We now look at the finite-size scaling functions for the spin-spin and energy-energy correlation 
lengths = N"^ S„^^{Ca^^z), see eq. (§). From universality with 5* = |, we expect 

27r 



Ls,(C.z) = Jl+'^ (12) 

27C V 77"^ 

where z = N{h — hc) is the finite-size scaling variable and -R^^i q(x) is a remnant function |T^. The 
spin-dependence should only enter into the metric factors Co- and which are determined from 
S„ and S^, respectively. In Figure 2, we display the extrapolated finite-size data of S„^^ for 77 = 0.7 
and find that they match nicely with the expected functional form. This confirms universality. 
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Table 3: Low lying excitations for charge Q = at the critical point. In each box, the upper value 
corresponds to 7] = 1, the lower one to r/ = 0.3. A dash indicates that no level is present, a '?' 
indicates that the finite-size data did not converge. For V — 0, all eigenstates shown have parity 
P — +1 and the lowest excitations with P — —1 occur for £ > 6. 
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Table 4: Low lying excitations for charge Q — 1. In each box, the upper value corresponds to 
T) — 1, the lower one to 77 = 0.3. 
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Table 5: Non- universal metric coefficients C^, found from the scaling fucntions and 5*^, 
respectively. 

Similar plots are obtained for other values of rj. The results for the metric factors are collected in 
table ||. Our results are consistent with 

CM = CM = Civ) (13) 

It is interesting to compare these with the conformal normalization ({rj) from table 0. Our data 
are roughly consistent with a linear relation C^^{ri) = a ■ (^{r]) with a ~ 0.75. 

A few comments are in order. Firstly, the observation eq. (|I3[) that the numerical value of 
the metric factor is independent of the physical quantity used for its determination, is certainly 
in agreement with the scaling expectation eq. (^) |1|]. Similar results were recently reported for 
2D percolation |jl6[|, where it was also checked that the metric factors are independent of the 
boundary conditions. Secondly, our results conffim earlier work [|1^ on the universality of the 
finite-size scaling function S„ in the 2D spin-1 Ising model. Thirdly, the observed linear relation 
between the conformal normalization ( and the metric factor C can be understood in terms of 
conformal perturbation theory, see ^ . In that framework, one would write for the non-critical 
quantum Hamiltonian H = ^{Hc + g4>), where He is the critical point quantum Hamiltonian, (p 
a perturbing relevant operator and g a non-universal coupling. In our case, (p = e, the energy 
density and g = h — h^. Since a given quantum Hamiltonian must in general be normalized to 
make conformal invariance applicable (see above), we note that into a perturbative calculation of 
the energy spectrum only the finite-size scaling variable N'^~^'i>{h — hc)/C enters. That is consistent 
with our finding C ~ 

In conclusion, we have investigated the ground state phase diagram of the spin-1 quantum XY 
chain in a transverse magnetic field. The structure of the phase diagram, obtained from finite-size 
scaling, is found to be very similar to the known spin-| case. We have explicitly conffimed the 
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universality of the Ising line, with respect to the spin S as well as the spin anisotropy rj, con- 
sidering both the conformal operator content and the finite-size scaling functions of the first two 
correlation lengths. 
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Figure captions 

Figure 1 Ground state phase diagram of the spin-1 XY model (|I|). P labels the disordered para- 
magnetic phase, F labels the ordered ferromagnetic phase and O labels the ordered oscillatory 
phase. The dotted line gives the P/F transition which for 7^ is in the 2D Ising universality 
class. The dashed line, given by 77 = 2yJ (h — l)/5 is the approximation to the P/F line as found 
from second order perturbation theory around rj = 0. The full line represents the F/0 transition 
as given by (|^). 

Figure 2 Finite-size scaling functions Sa{Caz) (lower curve) and Se{Cez) (upper curve) as a func- 
tion of the finite-size scaling variable z = N{h — he) for r/ = 0.7 as compared to the extrapolated 
finite-lattice estimates (points). 
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